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ABSTRACT 

The parsec-scale Faraday rotation measure (RM) distribution of six "blazars" 
(0954+658, 1156+295, 1418+546, 1749+096, 2007+777 & 2200+420) is investigated 
using multi-frequency (4.6, 5.1, 7.9, 8.9, 12.9, 15.4, 22 & 43 GHz) polarization obser- 
vations taken on 2006 July 2 with the American Very Long Baseline Array (VLB A). 
Analysis of the RM provides the direction of the line-of-sight (LoS) magnetic held 
component, as well as the intrinsic (unrotated) 2-D polarization distribution on the 
plane of the sky. Our results show that the magnitude of the core RM increases sys- 
tematically with frequency, and is well described by a power-law, where |RM core |oc v a . 
Our measured values of a vary from 0.9 to 3.8, providing information on the assumed 
power-law fall-off in the electron density (n e ) with distance from the central engine (r) 
for each source. RM gradients were detected across the jets of 0954+658, 1156+295 
and 1418+546, supporting the presence of helical magnetic fields in a sheath or bound- 
ary layer surrounding their jets. We find a bi-modal distribution of the intrinsic jet 
polarization orientation, with the polarization angles either aligned or orthogonal to 
the jet direction. The polarization of 2200+420 displays a continuous structure, with 
the polarization angles remaining aligned with the jet direction even as it bends. 
This indicates that the magnetic field structure in the synchrotron emitting plasma 
is dominated by an ordered transverse component. A helical magnetic field geometry 
can neatly explain both the bi-model distribution of the jet polarization orientation 
and the ordered polarization structure on these scales. For 0954+658, 1418+546 and 
2200+420, we find that the core RM changes sign with distance from the central en- 
gine. We provide an explanation for this by considering a boundary layer of Faraday 
rotating material threaded by a helical magnetic field, where bends in the relativis- 
tic jet or accelerating/decelerating flows give rise to changes in the dominant LoS 
components of the magnetic field, which in turn gives rise to different signs of the 
RM. 
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1 INTRODUCTION 



Active Galactic Nuclei (AGN) jets emit synchrotron radi- 
ation that can be detected at all radio frequencies. This 
radiation is often highly linearly polarized, which provides 
important information on the degree of order and orienta- 
tion of the magnetic field in these jets. The type of AGN 
we will consider here are known as "blazars", which have 
their jets pointed close to our line of sight (LoS) and often 
exhibit strong variability in total flux and linear polariza- 
tion over a broad range of frequencies from 7-ray to radio. 
Even though these AGN are intrinsically two-sided, pc-scale 
observations typically show a one sided "core-jet" structure 
due to Doppler boosting of the radiation from the relativis- 



tic jet pointed towards us. The jet moving away from us is 
usually too faint to be detected. 

Recent magneto-hydrodynamic (MHD) simulations 
have provided an almost complete theoretical picture of how 
these jets are launch ed, accelerate d and collimated close 
to the black hole, see iMeierl |2009l) and references therein. 
These models typically suppose a black hole that collects 
material in a magnetized, rotating accretion disk. A helical 
magnetic field is then produced by the differential rotation 
of the accretion disk or black hole ergosphere, and the mag- 
netocentrifugal forces produced by this field expel material 
out from the disk along the field lines, while magnetic pinch- 
ing forces collimate the jet outflow. It is possible that the 
magnetic field structure is randomised when the flow gets 
disrupted by shocks after a few hundred Schwarzschild radii, 
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but observ ational evidence of helical fields on scales larger 
than this llAsada et aT 2002; Gabuzda. Murray fc Croninl 
12004 : IZavala fc Tavlorl 120051 : flVIahmud fc Gabuzdal l2008h 
suggests that remnants of the earlier magnetic field struc- 
ture remain, or possibl y that a current driven helical kink 
instability is generated (Nakamura fc Meier 2004; Nak amural 
l2009l : ICarev. Sovinec fc Heinzll2008l ).~ 

The observed spectra of AGN jets can be optically thick 
(optical depth, r, greater than unity), self-absorbed (con- 
tribution from both optically thick and optically thin emis - 
sion regions), or optically thin (e.g.. lBeckert fc Fa lckc 2002). 
The radio cores of AGN are generally flat-spectrum, self- 
absorbed regions dominated by emission from around the 
t = 1 surface and the inner jet. The extended jet emission 
is generally optically thin, but can have regions of flatter 
spectra due to particle re-acceleration or interaction with 
the surrounding medium. Free-f ree absorption is sometimes 
detec ted at low frequencies (e.g.. lGabuzda. Gomez fc Agudol 
l200ll ). In the optically thin regime, the polarization orienta- 
tion is perpendicular to the magnetic field, with a theoretical 
maximum degree of polarization of ~ 75%; in the optically 
thick regime, the polarization is parallel to the magnetic 
field and th e degree of polari zation can be no more than 
~ 10 - 15% i|Pacholczvkll 1970h . 

Multi-frequency VLBI polarization observations allow 
the effect of Faraday rotation in the immediate vicinity of 
the AGN to be investigated. Faraday rotation of the plane 
of polarization of an electromagnetic wave occurs when it 
propagates through a region with free electrons and a mag- 
netic field. The effect is manifest as a linear dependence of 
the observed polarization angle (xobs) on the wavelength (A) 
squared, described (in SI units) by the formula 



Xobs = Xo + 



e 3 A 2 



nB ■ dl = xo + RMA 



(1) 



8n 2 EoTn 2 c 3 

l|jacksonlll975T ). where xo is the intrinsic (unrotated) polar- 
ization angle and the rotation measure, RM, is proportional 
to the integral of the Faraday rotating particle density (n) 
and the dot product of the magnetic field (B) with the path 
length (dl) along the LoS. Hence, a positive RM tells you 
that the LoS component of the magnetic field is pointing to- 
wards the observer and a negative RM means that the LoS 
component of the magnetic field is pointing away from the 
observer. 

For a pure electron-positron plasma, one would ex- 
pect zero net Faraday rotation because of the cancellation 
of the equal but opposite rotations from the two types of 
charged particle due to the e 3 dependence. The degree of 
Faraday rotation scales with the inverse square of the ef- 
fective mass (7m) -2 , so that the effect is much greater for 
thermal electrons tha n for relativistic electrons or protons 
l|jones fc Q'Delll977h . The observed Faraday rotation could 
be from the low energy end of the emitting electron popu- 
lation or from non-radiating thermal elect rons possibly en- 
trained in a sheath surrounding the jet ()Hughes fc Miller] 
1991). However, if the rotating medium is mixed in with the 
emitting plasma then internal Faraday rotation can occur. 
Theoretical models of inte rnal Faraday rotation in a spher - 
ical or cylindrical region (|Burnlll966l : ICioffi fc Jonel [l980) 
have shown that, if a rotation of greater than 45° is observed, 
then the Faraday rotation must be external . Many VLBI ob- 
servations show rotations greater than 45° (I Zavala fc Tavlorl 



I2003L |2004 iGabuzda et"al]|2004 120061 "). indicating that the 
main contribution to the rotation must be from gas segre- 
gated from the synchrotron e mitting region. Furtherm ore, 
observations of RM variability IjZavala fc Tavlorl 120011 ) and 
of transvers e RM gradients associated with helical m a gnetic 
fields (e.g., lAsada et all 120021 : IGabuzda etaD 12004 120081 ) 
indicate that the Faraday rotating gas must be very close 
to the source, most likely in a boundary layer or sheath 
surrounding the jet. Therefore, multi-frequency observations 
enable us to see gas that would otherwise be invisible to us 
and, combined with the corrected polarization orientation, 
which corresponds to the magnetic field projected onto the 
plane of the sky, we are provided with a three dimensional 
(3-D) view of the magnetic field structure. 

A number of previous studies of pc-scale Faraday rota- 
tion llTavlorl 1 19981. 120001: iRevnolds. Cawthorne fc Gabuzdal 
200 it IGabuzda fc Chernetskiil I2OO3I ; IZavala fc Tavlorl 120031 
200J) have focused on enhancements in magnitude of the 



RM in the VLBI core compared to values in the VLBI jet. 
Some previous results indicated the presence of different RM 
signs in t he core regions of "blazars" i n different frequency 
intervals (jO'Sullivan fc Gabuzdall2006h . Although RM vari- 
ability could not be ruled out in this case, the other possi- 
bility was that the direction of the dominant LoS magnetic 
field was changing with distance from the centre of activity. 
The 5-43 GHz VLBI observations, for which we present re- 
sults here, were designed to investigate the magnitude and 
sign of the RM in six AGN jets over a wide range of fre- 
quencies/spatial scales. 

Section 2 describes the observations and data reduction 
procedure. Our results are presented in Section 3, followed 
by their discussion in Section 4. Some conclusions are drawn 
in Section 5. We assume a cosmology with Ho = 71 km s _1 
Mpc" 1 , n M = 0.27 and Q A = 0.73. We define the spectral 
index, a, such that S y oc v +a . 



2 OBSERVATIONS AND DATA REDUCTION 

We present simultaneous multi-frequency polarization ob- 
servations obtained on the American Very Long Baseline 
Array (VLBA) radio interferometer, which provides milliarc- 
second (mas) resolution corresponding to the pc-scale struc- 
ture of these jets, over a 24-hour period on 2006 July 2. Six 
"blazars" (listed in Table [1} were observed at eight frequen- 
cies from 4.6 to 43 GHz (Table [2| in a "snapshot" mode, 
with roughly 8-10 scans of each object spread out over the 
time the source was visible with most or all of the VLBA an- 
tennas. This provided uniform u-v coverage for all sources. 
The total data rate for the observation was 128 Mbits s _1 . 
The initial calibration was performed with the radio astron- 
omy software package AIPS using standard techniques. Los 
Alamos was used as the reference antenna at all stages of 
the calibration. The instrumental polarizations ( "D-terms" ) 
were derived with the AIPS task LPCAL using CLEAN mod- 
els for individual features in the compact polarized source 
1156+295, which had good parallactic angle coverage over 
the observational period. After applying the solution to the 
data, plots of the real vs. imaginary cross-hand polarization 
data indicated that the "D-terms" had been successfully re- 
moved. 

The calibration of the electric vector position angles 
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Table 1. Observed sources 



Source (Other Name) Type Redshift Integrated RM 

(rad m — 2 ) 



0954 + 658 


BL Lac 


0.368 


-15 ±4 


1156 - 


h295 (4C +29.45) 


Quasar 


0.729 


-35 ± 1 


1418 - 


-546 


BL Lac 


0.152 


17 + 7 


1749 - 


1- 096 (4C +09.57) 


BL Lac 


0.320 


95 + 6 


2007- 


r 777 


BL Lac 


0.342 


-20 + 3 


2200 - 


h 420 (BL Lac) 


BL Lac 


0.069 


-205 ± 6 



Table 2. Observed frequencies with corresponding bandwidths 
and EVPA corrections 



Frequency Bandwidth Ax 
(GHz) (MHz) (deg) 



4.608 


16 7 


5.088 


16 11 


7.912 


32 90 


8.879 


32 84 


12.935 


32 62 


15.379 


32 89 


22.232 


32 32 


43.132 


32 88 



(EVPAs) was determined using integrated polarization mea- 
surements from the Very Large Array (VLA) polarization 
monitoring prograrrQ and the 26m radio telescope operated 
at the University of Michigan Radio Astronomy Observatory 
(UMRAO). The source 2200+420 was used for our polariza- 
tion position angle calibration. We used integrated polariza- 
tion measurements from the VLA on 2006 July 4 at 5, 8.5, 
22 and 43 GHz, and from UMRAO at 4.8 (July 3), 8.0 (June 

1 http:/ /www. aoc.nrao.edu/~smyers/calibration/ 
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Figure 1. Plot of integrated \ values, from VLA and UMRAO 
monitoring of 2200+420, vs. A 2 . We use 5.0 and 8.5 GHz VLA 
angles (squares) with 4.8, 8.0 and 14.5 GHz UMRAO angles (tri- 
angles) to find an integrated RM of —198 ± 4 rad m~ 2 . 



Table 3. VLA, UMRAO and VLBA data for the EVPA calibrator 
2200+420 



Telescope 


Frequency 


Date 


p Flux 


X 




(GHz) 




I'm Iv^i 


(deg) 


VLA 


5.0 


2006 July 7 


161 


-32.0 




8.5 


2006 July 7 


142 


-5.5 




22 


2006 July 7 


136 


10.5 




43 


2006 July 7 


106 


33.5 


UMRAO 


4.8 


2006 July 3 


165 


-36.0 




8.0 


2006 June 18 


143 


-13.5 




14.5 


2006 June 26 


132 


4.0 


VLBA 


4.6 


2006 July 2 


161 


-46.5 




5.1 


2006 July 2 


150 


42.5 




7.9 


2006 July 2 


128 


82.0 




8.9 


2006 July 2 


134 


-87.0 




12.9 


2006 July 2 


126 


-59.0 




15.4 


2006 July 2 


141 


-82.0 




22 


2006 July 2 


133 


-21.0 




43 


2006 July 2 


107 


-54.0 



18) and 14.5 GHz (June 26). These observations were com- 
pared with our total VLBI polarized flux (p = yj Q 2 + U 2 ) 
and position angle (x = |arctan^) for 2200+420 on 2006 
July 2. Since this calibration scheme relies on the fact that 
a large fraction of the polarized flux observed by the VLA 
is also detected with the VLBA (i.e., negligible intermediate 
scale polarized structure), it is very important for the po- 
larized fluxes to match. Table [3] shows that, for 2200+420, 
more than 90% of the VLA and UMRAO polarization was 
present on the VLBA scales at all frequencies. Corrections 
for the EVPAs were then found by aligning the VLA and 
VLBA polarization position angles (Ax = Xint — Xobs)- The 
polarization monitoring observations were interpolated for 
our closely spaced frequency pairs in accordance with the 
A 2 law displayed by the integrated EVPAs (see Fig. [TJ. The 
integrated RM indicated by these measurements is consis- 
tent with the v alue of —205 + 6 rad m~ 2 measured at longer 
wavelengths by iRudnick fc Jonesl (|l983f i. These corrections 
(Table [2J, relative to the same reference antenna, appear to 
be stable across long periods of time, and we estimate that 
they are accurate to within ±3°. 

The imaging was performed in AIPS using "robust" 
weighting, which is a combination of the natural and uniform 
weighting schemes. This weighting scheme is very useful for 
VLBA data of AGN jets because good resolution is needed 
to image the core region but we also want good sensitivity 
to the extended jet structure. Distributions of the polarized 
flux (p) and polarization angle (%), as well as accompany- 
ing "noise maps" , were constructed for all sources using the 
AIPS task COMB. The uncertainties written in the output 
X noise maps were calculated in COMB using the rms noise 
levels on the input Stokes Q and U maps. 

2.1 Analysis of the Rotation Measure 

In general, the total Faraday rotation for AGN is domi- 
nated by two components: (1) The fo reground contrib ution 
(mainly from our own Galaxy, see e.g.. |Pushkarevl (|200lh and 
references therein), and (2) magnetized thermal plasma in 
the immediate vicinity of the radio source. The foreground 
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Figure 2. Plots of \obs vs - ^ 2 f° r the whole observed frequency range of the core regions of the full resolution images of 0954+658 (a), 
1156+295 (b), 1418+546 (c), 1749+096 (d), 2007+777 (e) and 2200+420 (f). Dotted lines inserted where obvious linear dependences are 
present. Clear transitions in 0954+658, 1156+295, 1418+546 and 2200+420 outline where physically seperate core RM regions are being 
probed along the jet. The hollow data points for 2007+777 (e) are rotated by 90° due to the detection of an optical depth transition in 
the core between 8.9 and 12.9 GHz. The 4.6 and 5.1 GHz core EVPAs of 2200+420 (f) are not included in the plot in order to more 
clearly show the core RM transition at the shorter wavelengths; see Fig. 120b for the 4.6 and 5.1 GHz angles. 



© 2009 RAS, MNRAS 000, [Tp?? 



3-D Magnetic Field Structure in AGN Jets 5 



Table 4. Spectral indices corresponding to RM ranges 









C^corc 


Ojet 


"jot 


"jot 


Source 


Low v range 


Mid v range 


High v range 


Low v range 


Mid v range 


High v range 


0954 + 658 


0.3 


0.2 


0.2 


-0.5 


-0.6 


-0.5 


1156 + 295 


1.3 


0.5 


0.0 


-0.8 






1418 + 546 


0.4 


0.3 




-0.2 


-0.8 




1749 + 096 


0.8 


0.0 


-0.1 








2007 + 777 


0.6 


0.3 


0.2 


-0.4 


-0.6 


-0.6 


2200 + 420 


0.7 


0.5 


0.5 


-0.5 


-0.7 


-0.5 



Spectral indices for the low, mid and high v ranges arc derived between the highest and lowest frequencies in the corresponding RM ranges in 

Tablc[5] The errors in the spectral index values are ±0.1. 



Table 5. Summary of core RM results, with all RM values quoted in rad m 2 . 



Source 


RMcorc 


f^corc 


RMcorc 


W>corc 


RMcorc 




a 


Name 


(Low v range) 


(4.6 GHz) 


(Mid v range) 


(7.9 GHz) 


(High v range) 


(15.4 GHz) 


(| RM CO rc |0C V a ) 


0954 + 658 


+33 ± 14 


7.2% 


-88 ± 23 


5.1% 


-1591 ± 265 


1.2% 


3.84 + 1.34 


1156 + 295 


+170 + 5 


2.6% 


+618 ± 91 


2.8% 


+1667 ± 159 


1.3% 


2.22 + 0.05 


1418 + 546 


+83 ± 11 


4.3% 


-501 ± 48 


3.4% 




1.8% 


3.32 + 0.60 


1749 + 096 




3.3% 




2.9% 




3.3% 




2007 + 777 


+638 ± 39 


2.8% 




1.4% 


+1630 ± 201 


5.7% 


0.91 + 0.62 


2200 + 420 


-193 + 29 


1.3% 


+240 ± 90 


0.9% 


-1008 + 43 


3.6% 


1.40 + 0.18 



Low v range: generally 4.6—8.9 GHz; Mid v range: generally 7.9-15.4 GHz; High v range: generally 15.4-43 GHz. Exact ranges (sec text and RM 
maps) depend on where clear RM transitions occur and on the level of polarized flux detected in the jet at higher frequencies. No polarization 
was detected for 1418+546 at 22 and 43 GHz, hence there is no high f range RM for this source. No reliable Xobs vs. A 2 fits were found for 
1749+096 across the entire frequency range. m corc : core degree of polarization at the indicated frequency. The calculated power-law 
dependence, a, of the core RM with frequency for each source is quoted in the final column. 



contribution is constant over the entire radio source on milli- 
arcsecond scales, whereas the magneto-ionic Faraday rotat- 
ing medium near the radio source may be non-uniform on 
these scales; for example, the magnitude of the RM is usu- 
ally enhanced in the VLBI core region compared to its value 
in the VLBI jet (e.g.. IZavala fc Tavlodl2003ll2004l ). 

Removing the integrated RM is very important for our 
analysis for two main reasons. (1) The sign of the RM car- 
ries information about the direction of the LoS magnetic 
field giving rise to the Faraday rotation. In this sense, the 
sign of the RM contains information about the 3-D struc- 
ture of the magnetic field, which is lost if we are not able to 
isolate Faraday rotation local to the AGN from the Galactic 
component. (2) The relationship between the observed RM 
(RM bs) and the RM in the rest frame of the source (RMo) 
is RM obs = RM /(l + zf. For AGN with relatively low red- 
shifts, the values of RM b s and RMo will be very close, but 
for those with higher redshifts, it is essential to remove the 
foreground RM contribution if we are interested in compar- 
ing the intrinsic RMs in different AGN at different redshifts. 
Therefore, before making the RM maps with the AIPS task 
RM, we first removed the contribution of the known inte- 
grated Faraday rotation at each frequency, so that any re- 
maining RM should be solely due to magnetized thermal 
plasma in the vicinity of the AGN. 

Observations with both long and short frequency spac- 
ings are required for reliable Faraday rotation analysis, due 
to possible ±n7r ambiguities in the observed polarization an- 
gles. Our "broadband" frequency observations (5-43 GHz) 
enable us to construct RM maps on a range of different 
scales and Faraday depths, providing a wealth of informa- 
tion on the RM distribution and sign in the core region and 



out along the jet. It's important to note that VLBI resolu- 
tion is usually not sufficient to completely resolve the true 
optically thick core, therefore, the VLBI "core" consists of 
emission from the true core and some of the optically thin 
inner jet. Also, observations at different frequencies are de- 
tecting radiation near different r = 1 surfaces, so different 
frequency intervals probe different physical scales of the in- 
ner jet. Hence, matched-resolution images, convolved with 
the lowest frequency-interval beam, were constructed for dif- 
ferent frequency intervals where clear transitions in RM oc- 
curred (e.g., Fig. [2j . 

For the purposes of determining accurate RM values for 
individual regions in the VLBI core and jet, we found the 
mean polarization angle \ within a 3 x 3 pixel area at the 
same location in the polarization maps at each frequency. We 
identify the "core" as the furthest upstream self-absorbed jet 
region, which for some sources, especially at low frequencies, 
is upstream of the total intensity peak. We arbitrarily choose 
3x3 pixels (0.3 x 0.3 mas) as the area used to obtain the 
EVPA value and an estimate of its rms deviation in the cor- 
responding x noise map, which is then added in quadrature 
with the estimated EVPA calibration uncertainty of 3°. Fi- 
nally, a reduced % 2 linear regression was performed to find 
the RM and its standard error. 



3 RESULTS 

We first outline some general methods and results before 
considering each individual source. 
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3.1 Spectral Index 

Spectral index maps (S v on v a ) were constructed between 
pairs of frequencies to obtain detailed spectral information 
across the entire frequency range. Each map was corrected 
for the frequency dependence of the optically thick core po- 
sition using a cross-correlation tec hnique for aligning the op - 
tically thin emission from the jet l|Croke & Gabuzda 2008). 
This method is extremely useful for smooth or weak jets 
that do not have distinct components at each frequency 
with which to align the images. For sources with very 
sparse or negligible jet emission, in particular 1156+295 and 
1749+096 at high frequencies, we are unable to use this 
alignment method. In the case of 1156+295, where we have 
calculated core-shifts at lower frequencies, we predicted the 
shifts expected at higher freque ncies assuming that equipar- 
tition holds in these regions (see[Lpbanov 1993) an d checked 
that the resulting spectral index maps were consistent. For 
1749+096, we simply aligned the images by their peak flux 
densities and checked that the spectral index maps seemed 
reasonable. 

All sources display a flat or inverted spectrum across 
all frequencies in the furthest upstream region, consistent 
with a self-absorbed core. In general, the jet emission for 
all sources was optically thin (a ^ —0.5), with some 
flat spectrum regions possibly associated with particle re- 
acceleration or interaction with the surrounding medium. 
Detailed information on the spectra is required in order to 
identify any regions where possible optically thick-thin tran- 
sitions occur, requiring the observed polarization angles to 
be rotated by 90° for correct Faraday rotation analysis. A 
detailed analysis of our spectral index maps will be presented 
in a forthcoming paper. Spectral indices for both core and 
jet regions corresponding to the RM ranges for each source 
are summarized in Table [3] 



3.2 Rotation Measure 

In the current paradigm, Faraday rotation in the vicinity of 
the AGN occurs in a boundary layer or sheath of magnetised 
plasma surrounding the jet (Zavala & Taylor 2005, Attridge 
et al. 2005). Our observed RMs are consistent with the ro- 
tating medium being external to the emitting region. Linear 
fits of x vs - ^ 2 are found in most cases, with rotations of 
greater than 45° observed in several cases. Furthermore, the 
distribution of the degree of polarization of those with rota- 
tions less than 45° does not fit the corresponding predicted 

depolarization expected for intern al Faraday rota tion. 

Recent observations dGabuzda et all l2006t 
lO'Sullivan fc Gabuzda! 120061 ; Ijorstadl I2007T I have sug- 
gested that the magnitude of the observed core RMs in 
AGN increase systematically toward higher frequencies. 
Since higher frequency observations probe regions closer to 
the base of the jet, higher magnitude RMs simply imply that 
the electron density and/or the LoS magnetic field strength 
is increasing closer to the central engine. The frequency de- 
pendence of the core RM can be inve stigated directly using 
the relation derived in lJorstadl (|2007l ). |RM coro |cx v a , where 
a describes an assumed power-law fall-off in the electron 
density with distance r from the central engine, n e oc r~ a . 
We assume that the fitted core RM corresponds to the core 
region defined in the lowest frequency observation used 



(e.g., |RM C ore,!/ 1 -!/ 2 | tx (^2) a ,f2 < Vi). A summary of our 
core RM results is presented in Table [5] along with our 
estimated values of a. 



3.3 Intrinsic Polarization Orientation 

After correction for Faraday rotation, we are able to find the 
intrinsic polarization orientation of the core and jet; hence, 
we can make firm deductions on the intrinsic magnetic field 
structure. It's important to note that the intrinsic polar- 
ization orientations can vary with frequency-interval simply 
because different RMs may imply different zero-wavelength 
(Xo) angles. To determine the orientation of the furthest up- 
stream polarized component that corresponds to the radio 
core, we use either the highest frequency-interval xo value 
or the 43-GHz core polarization angle. 



3.4 Results for Individual Sources 

We will now consider the results for each source in turn. In 
each case, maps are presented for six of our eight frequen- 
cies; the 5.1 and 8.9 GHz images are omitted due to their 
similarity to the 4.6 and 7.9 GHz images, respectively. 



0954 +658 

This source is classified as a BL Lac object with a redshift 
of 0.3 67. It's a well known Intra-Day Variable (IDV) source 
(e.g., iGabuzda et~ai1 |200Ch and has also been detected at 
•7-ravs by EGRET l|Mattox et al1l200ll ). It displays a one- 
sided kpc-scale jet structure extending to the southwest 
l|Cassaro et al.l Il999l ) , which seems consistent with a con- 
tinuation of the observed bend on pc-scales in the 4.6-GHz 
total intensity map (e .g., Fig. [3^,) . It has an integrated RM 



of -15 ± 4 rad m~' z (|RuskHl988T ). We find the VLBI core 
to be self-absorbed with a ~ +0.3 from 4.6-15.4 GHz and 
a ~ +0.2 from 15.4-43 GHz. The jet emission is generally 
optically thin with a ~ —0.5 from 4.6-15.4 GHz and from 
15.4-43 GHz (Table |4]|. A notable exception is the detec- 
tion of a gradient in the 5.1-7.9 GHz spectral index map 
(Fig. 2K) where a changes from negative to positive across 
the jet close to where it bends, possibly indicating interac- 
tion with the surrounding medium. Figure [3] displays the 
total intensity and polarization structure of 0954+658. 

The core degree of polarization (m) at 4.6 GHz is ~ 7%; 
it then decreases to a minimum of ~ 1.2% at 15.4 GHz before 
increasing again back up to ~ 7% at 43 GHz. The degree of 
polarization in the jet is higher and also has a clear gradient 
across the jet going from ~ 10% at the southern edge to 
~ 30% at the northern edge (Fig. Hp). This gradient in m 
is also seen on scales of 10 -20 mas in the same direction 
(|Hallahan fc Gabuzdj|20Qgt ). 

For this source, we were able to produce RM maps on 
three different scales: from 4.6-15.4, 7.9-15.4 and 12.9-43 
GHz. The RM map from 4.6-15.4 GHz (Fig. Hi) has a core 
RM of 33 ± 14 rad m" 2 . There is a clear gradient in RM 
across the inner jet region before the jet bends; the sign 
of the RM remains positive across the gradient, varying 



systematically from ~ 70 rad m 



to 



140 rad m"^ at 



the northern edge. The magnitude of the RM in the cores 
of AGN is usually larger than what is found in their jets 
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Figure 3. Total intensity (I) and polarization (p) maps for 0954+658: (a) 4.6 GHz; / peak 0.77 Jy bin -1 ; lowest contour 1.31 mjy 
bm" 1 ; p vectors clipped at 1.1 mjy bm" 1 , (b) 7.9 GHz; I peak 0.92 Jy bm -1 ; lowest contour 1.84 mjy bm~ 1 ; p vectors clipped at 1.1 
mjy bm -1 , (c) 12.9 GHz; / peak 1.02 Jy bm -1 ; lowest contour 2.24 mjy bm -1 ; p vectors clipped at 1.3 mjy bm -1 , (d) 15.4 GHz; I 
peak 0.95 Jy bm -1 ; lowest contour 2.38 mjy bm -1 ; p vectors clipped at 1.3 mjy bm -1 , (e) 22 GHz; / peak 0.93 Jy bm -1 ; bottom 
contour 5.37 mjy bm -1 ; p vectors clipped at 2.5 mjy bm -1 , (f) 43 GHz; I peak 1.04 Jy bm -1 ; lowest contour 8.81 mjy bm -1 ; p vectors 
clipped at 3.5 mjy bm -1 . The / contours increase by factors of two in all maps. Length of p vectors represent relative polarized intensity. 
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Figure 4. (a) Spectral index map for 0954+658 from 5.1—7.9 GHz (o range: -1.6 to 1.1), overlaid with 5.1 GHz total intensity contours 
(/ peak 0.78 Jy bm _1 ; lowest contour 1.63 mjy bm" 1 ). The slice displays the gradient across the jet that may indicate interaction with 
the surrounding medium where the jet bends, (b) 4.6 GHz total intensity (/ peak 0.77 Jy bm -1 ; lowest contour 1.61 mjy bm -1 ) and 
degree of polarization (m), ranging from to 50 percent, for 0954+658 showing how m increases near the region where the jet bends. 
The / contours increase by factors of two in both maps. Full resolution image in journal version. 



I Zavala fc Tavloij 120031 120041 ). but in this case the magni- 
tude is lowest in the core. This supports the above evidence 
from the spectral index map of possible interaction with the 
surrounding medium where the jet bends, with the higher 
electron density in this region causing the observed increase 
in the magnitude of the RM and the increase of the de- 
gree of polarization possibly due t o ordering of th e mag- 
netic field at the site of interaction (jAlov et al.ll200oh . How- 
ever, the RM and degree of polarization gradients can also 
be explained in term s of a helical magnetic field structure 
ijLvutikov et a.l.ll2005l ). This is an example of the difficulty 
to distinguish between asymmetric RM gradients associated 
with helcial magnetic fields and/or i nteraction with the sur- 
rounding medium in some cases. See[G omcz et alJ (|2008T I for 
a similar scenario in the radio galaxy 3C 120. 

Analyzing the RM from 7.9-15.4 GHz (Fig. Eh), we find 
that the core RM changes sign and increases in magnitude 
to -88 + 23 rad m~ 2 . The inner jet shows some evidence for 
a continuation of the RM gradient found at lower resolution, 
however the errors in the fits are very large in this region. 

The highest frequency-interval RM map, from 15.4-43 
GHz (Fig.[5t), only contains information about the core RM, 
since no appreciable polarization was detected outside the 
core region at 22 and 43 GHz. The core RM of -1591 ± 265 
rad m -2 has the same sign as the 7.9-15.4 GHz value but 
its magnitude has increased dramatically. Fitting the three 
core RM values vs. frequency using a least-squares method, 
we find a large value of a = 3.84 ± 1.34 for the fall-off in 
electron density from the central engine (n e oc r~ a ). 

After correcting for Faraday rotation, the intrinsic po- 
larization orientation derived from 4.6-15.4 and 7.9-15.4 



GHz, in both the core and jet, is approximately aligned 
with the jet direction (Fig. [6^., b). However, at the high- 
est frequency interval, 15.4-43 GHz, the furthest upstream 
polarized component is approximately transverse to the jet 
direction (Fig. [6]:). The 43-GHz polarization map clearly re- 
solves this component (Fig. [3f ) . 



1156+295 

Also known as 4C +29.45, this object has been classified 
as both a High Polarization Radio Q uasar (HPR Q) and an 
Optical Violent Variable (OVV) fsee iHond l20oj ) . It is the 
most distant object in our sample of six sources, with a 
redshift of 0.729. It has also displayed intra-day variability 
|Savolainen fc Kovalevl2008l) and has a 7-ray detection from 
EGRET (|Mattox et alj|200ll). On kpc-scales it has a two- 
sided north-south halo ()Antonucci fc Ulvestadl[l985l ) ; on pc- 
scales the jet initially propagates directly north and then 
after a few mas changes to a northeasterly direction, as can 
be seen in our 4.6-GHz image (Fig. [7 ^)- We subtracted an 
integrated RM of -35 ± 1 rad m -2 (|Oren fc Wolfd Il995l ) 
before making the RM maps. From our data, the spectral 
index in the core is inverted at low frequencies with a ~ 
+ 1.3 from 4.6-8.9 GHz, indicating strong synchrotron self- 
absorption. The inner jet is optically thin with a ~ —0.8. 
The core spectrum flattens at higher frequencies with a ~ 
0.0 between 12.9 and 43 GHz (Table |4j. Figure [7] displays 
the total intensity and polarization structure of 1156+295. 

X vs. A 2 fits for the core region were found in three sep- 
arate frequency ranges: 4.6-8.9, 7.9-15.4 and 12.9-43 GHz 
(Fig. [8]l . The core RM increases systematically with fre- 
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Figure 5. Source: 0954+658. (a) 4.6-15.4 GHz RM map (colour scale RM range: -20-170 rad m~ 2 ) with 4.6 GHz total intensity contours 
(/ peak 0.77 Jy bm _1 ; lowest contour 1.61 mjy bm _1 ). The solid arrow points to the core region where the \ vs - ^ 2 fits shown in the 
inset were found; the degree of polarization is plotted underneath for each corresponding polarization angle. The dashed arrow points 
to the indicated slice taken across the RM map; the plotted black line shows the RM fit at every pixel across the slice with the grey 
area indicating the error in the fit for each point, (b) 7.9—15.4 GHz RM map (colour scale RM range: -150—200 rad m -2 ) with 7.9 GHz 
total intensity contours (/ peak 0.90 Jy bm" 1 ; bottom contour 2.53 mjy bm" 1 ). (c) 12.9-43 GHz RM map (colour scale RM range: 
-3-0.5xl0 3 rad m~ 2 ) with 12.9 GHz total intensity contours (/ peak 0.98 Jy bm^ 1 ; lowest contour 2.76 mjy bm -1 ). All / contours 
increase by factors of two. Full resolution image in journal version. 
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Figure 6. Total intensity (I) contours with intrinsic polarization 
orientation (xo) f° r 0954+658 predicted from (a) 4.6—15.4 GHz 
RM map, (b) 7.9-15.4 GHz RM map and (c) 12.9-43 GHz RM 
map. Same / contour levels as in Fig. IB1 



quency from 170 ± 5 rad m~ 2 at 4.6-8.9 GHz, to 618 ± 91 
rad m~ 2 at 7.9-15.4 GHz and finally to 1667 ± 159 rad m -2 
at 12.9-43 GHz. This implies a fall-off in the electron den- 
sity from the central engine of n e oc r - 2 - 22 ± 05 . There is 
no detection of a change in sign of the core RM in different 
frequency intervals. The only jet RM, derived at the lowest 
frequency interval, of 136 ± 4 rad m -2 , is smaller in magni- 
tude than the RM in the core. There is evidence for a RM 
gradient across the inner jet in all three RM maps, as shown 
through the plots of the sliced regi ons in Fig. [8] confirming 
the results of iGabuzda et ail (|200Sf ). The gradient increases 
from ~ 50 rad m -2 at the eastern edge to ~ 250 rad m~ 2 
at the western edge of the jet in the 4.6-8.9 GHz RM map, 
from ~ 250 rad m~ 2 to ~ 650 rad m -2 in the 7.9-15.4 GHz 
map and from ~ 1300 rad m~ 2 to ~ 1800 rad m -2 in the 
12.9-43 GHz map. 

The core degree of polarization decreases from ~ 2.8% 
at 43 GHz to a minimum of ~ 1% at 15.4 GHz and back up 
to ~ 3% at 5.1 GHz with a slight decrease to ~ 2.6% at 4.6 
GHz. The degree of polarization is ~ 15% in the inner jet 
at 5.1 and 4.6 GHz. 

The outer jet polarization in the 4.6-GHz map displays 
a spine-sheath type structure with polarization aligned with 
the jet direction in the central region and transverse to it 
along the edges of the jet (Fig. [7^). In the inner jet region, 
the Faraday corrected polarization vectors are transverse to 
the jet direction (Fig. [UJi) but this is likely to correspond to 
the sheath polarization structure. The intrinsic core EVPA 
derived from the highest frequency interval (12.9-43 GHz) 
is ~ —20° which appears to be slightly offset from the jet 
direction (Fig. [9J;). However, a jet direction of ~ —20° at 
higher frequencies would n ot be inconsi stent with the helical 
jet trajectory proposed by iHon el l|2004h . 

1418+546 

This BL Lac object, at a redshift of 0.152, has a one-sided 
kpc-scale halo with a component west of the compact core 
l|Cassaro et al.|[l999h and is misaligned with the VLBI pc- 
scale jet, which propagate s in a south-easterly direction 
(Fig. [TO)). iRudnick fc Jones] (|l983h found an integrated RM 
of 17 ± 7 rad m~ 2 for this source. We find a core spec- 
tral index of a ~ +0.4 consistent with self-absorption and 
a jet spectral index of a ~ —0.2 from 4.6-8.9 GHz. The 
core region is resolved into two components at higher fre- 
quencies and the spectral index map from 7.9-15.4 GHz 
gives a ~ +0.3 for the furthest upstream component, with 
a ~ —0.6 and a ~ —0.8 for the inner and outer jet regions, 
respectively (Table [3). Figure [TOl displays the total intensity 
and polarization structure of 1418+546 at four of the six 
frequencies at which polarization was detected. 

We could only obtain RM maps for two frequency in- 
tervals for this source (4.6-8.9 and 7.9-15.4 GHz) due to a 
failure to detect polarization at 22 and 43 GHz. The core RM 
from 4.6-8.9 GHz is 83 ± 11 rad m -2 ; a RM of 91 ± 13 rad 
m -2 was found in the jet from the same frequency range 
(Fig. I lib ,') . There is also a gradient in RM across the jet 
from ~ —30 rad m -2 at the north eastern edge to ~ 120 
rad m~ 2 at the south western edge fFig. lllfa '). In the higher 
frequency interval, from 7.9-15.4 GHz, the RM in the core 
changes sign and increases in magnitude to —501 + 48 rad 
m -2 , giving an estimate of a ~ 3.3. The RM in the inner 
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Figure 7. Total intensity (/) and polarization (p) maps for 1156+295: (a) 4.6 GHz; / peak 1.05 Jy bin -1 ; lowest contour 2.95 mjy 
bm" 1 ; p vectors clipped at 1.1 mjy bm" 1 , (b) 7.9 GHz; I peak 2.16 Jy bm -1 ; lowest contour 4.55 mjy bm -1 ; p vectors clipped at 2.5 
mjy bm -1 , (c) 12.9 GHz; / peak 2.98 Jy bm -1 ; lowest contour 8.94 mjy bm -1 ; p vectors clipped at 2.8 mjy bm -1 , (d) 15.4 GHz; I 
peak 2.91 Jy bm -1 ; lowest contour 8.72 mjy bm -1 ; p vectors clipped at 2.8 mjy bm -1 , (e) 22 GHz; J peak 3.09 Jy bm" 1 ; bottom 
contour 17.62 mjy bm -1 ; p vectors clipped at 3.5 mjy bm" 1 , (f) 43 GHz; I peak 2.85 Jy bm -1 ; lowest contour 24.18 mjy bm -1 ; p 
vectors clipped at 5.5 mjy bm -1 . The I contours increase by factors of two in all maps. Length of p vectors represent relative polarized 
intensity. 



jet has a value of 61 ± 36 rad m -2 and a zero RM (within 
the errors) is found in the outer jet region (Fig. Illb ); there 
seems to be an apparent RM gradient in this map, however 
the errors in the higher frequency angles are large and the 
gradient is not consistently present along the outer jet, un- 
like the lower frequency RM gradient in approximately this 
same location. Hence, we don't consider this to represent a 
real RM gradient, only patchiness in the RM distribution. 

The degree of polarization in the core increases steadily 
from ~ 1.8% at 15.4 GHz to ~ 4.3% at 4.6 GHz. There is a 
gradient in m across the jet ranging from ~ 5% at the south 
western edge to ~ 25% at the north eastern edge (Fig.[T2"R). 
Since the magnitude of the RM gradient increases in the op- 
posite direction across the jet, and the spectral index map 



shows no transverse structure in this region (Fig. 112b ). in- 
teraction with the surrounding medium is unlikely to be a 
dominant effect in this case. This provides strong support 
for the direct association of the transverse RM gradient with 
the underlying magnetic field geometry. 

The dominant Faraday corrected intrinsic jet EVPAs 
are transverse to the jet direction (Fig. 113b . b) . The polar- 
ized emission in the outer jet is clearly offset to the north- 
eastern edge of the jet from the total intensity emission (e.g., 
Fig. I10[) . This is similar to the polarizati o n stru cture in the 
5-GHz observations of IPushkarev et al.l (|2004h . where the 
jet polarization is also offset to the northern edge of the jet. 
(However, note that, at another epoch, they observed the 
same jet region to have polarization offset to the southern 
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Figure 8. Source: 1156+295. (a) 4.6-8.9 GHz RM map (colour scale RM range: -70-320 rad m~ 2 ) with 4.6 GHz total intensity contours 
(I peak 1.05 Jy bm -1 ; lowest contour 2.73 mjy bm _1 ). The solid arrow points to the core region where the \ vs - ^ 2 shown in the 
inset were found; the degree of polarization is plotted underneath for each corresponding polarization angle. The dashed arrow points to 
the indicated slice taken across the RM map; the plotted black line shows the RM fit at every pixel across the slice with the grey area 
indicating the error in the fit for each point, (b) 7.9-15.4 GHz RM map (colour scale RM range : 0-1 xlO 3 rad m~ 2 ) with 7.9 GHz total 
intensity contours (I peak 2.16 Jy bm -1 ; lowest contour 4.33 mjy bm -1 ). (c) 12.9-43 GHz RM map (colour scale RM range: 0— 2.5xl0 3 
rad m — 2 ) with 12.9 GHz total intensity contours (/ peak 2.91 Jy bm -1 ; lowest contour 7.27 mjy bm -1 ). All / contours increase by 
factors of two. Full resolution image in journal version. 
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Figure 9. Total intensity (/) contours with intrinsic polarization orientation (xo) f° r 1156+295 predicted from (a) 4.6—8.9 GHz RM 
map, (b) 7.9-15.4 GHz RM map and (c) 12.9-43 GHz RM map. Same / contour levels as in Fig. [8] 



edge of the jet at both 5 and 8.4 GHz.) From 4.6-8.9 GHz, 
the intrinsic EVPAs in the core are also approximately trans- 
verse to the jet direction (Fig. 113a ). However, the higher fre- 
quency interval (7.9-15.4 GHz) core xo shows polarization 
in the furthest upstream component approximately aligned 
with the jet direction (Fig. 113b ). 



1749+096 

This BL Lac object has a redshift of 0.32 2 and is a point 
source on kpc scales l|Rector fc Stockel200ll ). On pc scales, it 
has a jet that extends in a north-easterly direction (Fig. 1 141) . 
althou gh the VSOP polarization observations of iGabuzdal 
(|2005f ) show a new jet component emerging in a structural 
position angle of ~ —15°. This s ource has a relat ively large 
integrated RM of 95±6 rad m" 2 |Pushkarevl200ll ). The core 
region shows clear evidence for a spectral turnover, with the 
spectral index decreasing from a ~ +0.8 between 4.6-8.9 
GHz toa~ -0.1 between 12.9 and 43 GHz (TablegJ. Figure 
[redisplays the total intensity and polarization structure of 
1749+096. 

Reduced \ 2 fits of the observed core polarization an- 
gles vs. A 2 produced no acceptable fits across our entire 
frequency range for this source despite several attempts 
involving ±7i7r rotations and ±90° rotations where signif- 
icant changes in the spectral index were observed (the lat- 
ter due to possible optically thick-thin transitions); see 
Fig. [15] for a plot of Xobs vs. A 2 . No clear evidence of a 90° 
jump between the polarization angles at adjacent frequen- 
cies is observed, despite the evidence for a spectral turnover 
from the measured core spectral indices. This is surpris- 
ing since previous observations have yielded good RM fits 



for this source (IZavala fc Tavlonbooi iGabuzda et~aL1l2006l : 



iMahmud fc Gabuzdal 20081 ). We may have observed the 
source during a flare since the peak flux for our observations 
is much larger than what has previously been observed. For 
example, our peak flux density at 7.9 GHz was 4.2 Jy/beam 
(epoch 2006), with previous observations that did detect 
Faraday rotation (epoch 2004) having a peak flux density 
of 3.0 Jy/beam at 7.9 GHz (Mahmud, priv ate communica- 
tion) ; also, previous 43-GHz observations (Gabuzda et al. 
2006) have a peak flux density of 2.4 Jy/beam whereas 
the peak flux density in our 43-GHz m ap is 3.9 Jy/beam. 
ISokolov. Cawthorne fc O'Sullivanl (|2008l ) detected a corre- 
lation between the polarized flux in the 43-GHz core and 
the 15-43 GHz core RM in multi-epoch observations of BL 
Lac indicating significant interaction between the emitting 
plasma and the Faraday rotating material. Hence, if the ro- 
tating material is mainly in a boundary layer and if it was 
affected by a flare, this may explain the lack of a clear A 2 fit 
for 1749+096 at this epoch. 



The core degree of polarization remains relatively con- 
stant across the whole frequency range, with m ~ 3% 
(Fig. 1151 bottom) consistent with the birth of a new polar- 
ized component that dominates the emission. The absence 
of depolarization also rules out the effect of internal Fara- 
day rotation in the core being responsible for the lack of a 
A 2 dependence with polarization angle. There was no ap- 
preciable jet polarization detected in this source. The core 
polarization orientation is approximately transverse to the 
visible VLBI jet at all frequencies (Fig. I14|l . 
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Figure 10. Total intensity (/) and polarization (p) maps for 1418+546: (a) 4.6 GHz; I peak 0.56 Jy bm" 1 ; lowest contour 1.40 mjy 
bm _1 ; p vectors clipped at 1.0 mjy bm" 1 , (b) 7.9 GHz; / peak 0.59 Jy bm -1 ; lowest contour 1.38 mjy bm -1 ; p vectors clipped at 1.0 
mjy bm -1 , (c) 12.9 GHz; / peak 0.53 Jy bm -1 ; lowest contour 1.75 mjy bm -1 ; p vectors clipped at 1.5 mjy bm -1 , (d) 15.4 GHz; I 
peak 0.39 Jy bm -1 ; lowest contour 1.96 mjy bm -1 ; p vectors clipped at 1.5 mjy bm -1 . The / contours by factors of two in all maps. 
Length of p vectors represent relative polarized intensity. 



2007+777 

This source is classified as a BL Lac objec t with a redshift 
of 0.3 42. It has a two-sided kpc-scale jet (|Antonucci et all 
1986) orientated in the same direction as the one-sided VLBI 
pc-sca le jet (Fig. I16[) . The integrated RM found by Rusk 
(1988) is -20 ± 3 rad m~ 2 . Figure [16] displays the total 
intensity and polarization structure of 2007+777. The core 
is self-absorbed with a ~ +0.6 from 4.6-8.9 GHz, while the 
jet is optically thin with a ~ -0.4. From 12.9-43 GHz the 
core spectral index decreases to a ~ +0.2 and a ~ —0.6 
in the jet (Table 0}. The plot of Xobs vs. A 2 in Fig. [2] is 
suggestive of an optically thick-thin transition between 8.9 
and 12.9 GHz. This is supported by the large increase in the 
core degree of polarization between these frequencies (see 
Table [5}, as well as the flattening of the core spectrum from 



a ~ 0.6 between 4.6 and 8.9 GHz to a ~ 0.2 from 12.9- 
43 GHz. Therefore, we should rotate the core EVPAs from 
4.6 to 8.9 GHz by —90° before comparing with the higher 
frequency EVPAs (see Fig. [2^) . 

We constructed RM maps from 4.6 to 8.9 GHz and from 
12.9 to 22 GHz. The 43-GHz core EVPA does not lie on the 
A 2 line from 12.9 to 22 GHz, indicating that a physically 
different region is being probed. From 4.6-8.9 GHz, the RM 
map has a sharp boundary between the core and the inner 
jet (Fig. 117b ); this sharp boundary is also present in the 
spectral index map. The core has an RM of 638 + 39 rad m~ 2 
while the inner jet RM is lower at 310 + 41 rad m -2 . Out in 
the jet the RM is lower again at 99 ± 26 rad m -2 . From the 
12.9-22 GHz frequency interval, the core RM increases to 
1630 + 201 rad m -2 giving an estimate of a ~ 0.9. The inner 
jet has a RM of ~ rad m -2 within the errors (Fig. 117b ). 
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Figure 11. Source: 1418+546. (a) 4.6-8.9 GHz RM map (colour scale RM range: -100-160 rad m -2 ) with 4.6 GHz total intensity 
contours (I peak 0.55 Jy bm" 1 ; lowest contour 1.66 mjy bm _1 ). The solid arrow points to the core region where the x vs - fits shown 
in the inset were found; the degree of polarization is plotted underneath for each corresponding polarization angle. The dashed arrow 
points to the indicated slice taken across the RM map; the plotted black line shows the RM fit at every pixel across the slice with the 
grey area indicating the error in the fit for each point, (b) 7.9-15.4 GHz RM map (colour scale RM range: -1-0.4X10 3 rad m — 2 ) with 
7.9 GHz total intensity contours (I peak 0.59 Jy bm -1 ; lowest contour 1.76 mjy bm -1 ). All / contours increase by factors of two. Full 
resolution image in journal version. 
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Figure 12. (a) Map of the degree of polarization (m), ranging from to 30 percent, at 4.6 GHz for 1418+546 showing a gradient across 
the jet (I peak 0.56 Jy bm" 1 ; lowest contour 1.68 mjy bm" 1 ). (b) Map showing spectral index (o range: -1.5 to 1) between 5.1 and 7.9 
GHz, with 5.1 GHz contours (I peak 0.56 Jy bm -1 ; lowest contour 1.69 mjy bm -1 ). 
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Figure 13. Total intensity (/) contours with intrinsic polarization orientation (xo) for 1418+546 predicted from (a) 4.6-8.9 GHz RM 
map and (b) 7.9—15.4 GHz RM map. Same / contour levels as in Fig. 1111 



The core degree of polarization decreases rapidly from 
~ 6.5% at 43 GHz to ~ 0.9% at 8.9 GHz before rising slowly 
to ~ 2% at 4.6 GHz. Where jet polarization is detected at 
lower frequencies, the degree of polarization remains rela- 
tively constant (m ~ 8-10%). 



The Faraday corrected jet polarization is aligned with 
the jet direction in the 4.6-8.9 GHz xo map, while the core 
Xo does not reflect the true xo distribution in this region 
because the polarized emission is dominated by optically 
thick regions at these frequencies (Fig. IPHk). The inner jet 
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1749+096 




Figure 14. Total intensity (/) and polarization (p) maps for 1749+096: (a) 4.6 GHz; / peak 2.61 Jy bm" 1 ; lowest contour 2.09 mjy 
bm _1 ; p vectors clipped at 2.5 mjy bm -1 , (b) 7.9 GHz; / peak 4.19 Jy bm -1 ; lowest contour 3.77 mjy bm -1 ; p vectors clipped at 2.5 
mjy bm -1 , (c) 12.9 GHz; / peak 4.54 Jy bm -1 ; lowest contour 4.08 mjy bm -1 ; p vectors clipped at 3.0 mjy bm -1 , (d) 15.4 GHz; I 
peak 3.97 Jy bm -1 ; lowest contour 5.56 mjy bm -1 ; p vectors clipped at 3.0 mjy bm _I , (e) 22 GHz; J peak 3.87 Jy bm -1 ; bottom 
contour 9.67 mjy bm -1 ; p vectors clipped at 4.5 mjy bm -1 , (f) 43 GHz; I peak 3.86 Jy bm -1 ; lowest contour 38.62 mjy bm -1 ; p 
vectors clipped at 5.0 mjy bm" . The / contours increase by factors of two in all maps. Length of p vectors represent relative polarized 
intensity. 



and core xo values implied from the 12.9-22 GHz RM map 
are aligned with the jet direction (Fig. 118b ). The 43-GHz 
core EVPA does not follow the A 2 dependence from 12.9- 
22 GHz, possibly indicating that a physically separate RM 
region is being probed at this frequency. 



2200+420 

This widely studied BL Lac object (see Bach 2006, and ref- 
erences therein) is the nearest source in our sample, with a 
redshift of 0.0686. It has been detected in 7-rays by EGRET 



jMattox et alj 200 if) a nd it has a kpc-scale core-halo struc- 
ture <|Antonuccilll986h . 



It has a large inte grated RM of —205 ± 6 rad m 2 
|Rudnick fc Jonea[l9 83). The core region is self-absorbed, 
with a ~ +0.7 from 4.6-12.9 GHz and a ~ +0.5 from 15.4- 
43 GHz. The extended jet emission is optically thin, with 
q ~ —0.5 (Table 2J) . Figure [191 displays the total intensity 
and polarization structure of 2200+420. 

Three separate RM maps were constructed from 4.6- 
12.9, 7.9-12.9 and 15.4-43 GHz. There was a clear transi- 
tion in the A 2 dependence of the core EVPA between 12.9 
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Figure 15. Plot of Xobs vs. A 2 for core of 1749+096. There is no 
evidence for a linear dependence indicating external Faraday ro- 
tation despite several attempts involving nir rotations and ±90° 
rotations where significant changes in the spectral index were ob- 
served. The plot of the degree of polarization (m) underneath 
rules out the effect of internal Faraday rotation since m remains 
approximately constant across the whole frequency range. 



and 15.4 GHz, indicating a physically separate region of the 
jet was being probed (e.g., Fig. [2j bottom right). Interest- 
in gly, this transition also a ppears in the RM maps of BL Lac 
bv lZavala fc Taylor! (|2003l . Fig. 25) but without the 22 and 
43-GHz data it could not be identified. After removing the 
integrated RM, we get a value of 39 ± 9 rad m -2 out in the 
jet from the low resolution (4.6-12.9 GHz) RM map. This 
is consistent with a low electron density, as expected for the 
extended emission. The co re RM of —193 ± 29 r ad m~ 2 is 
similar to previous values (|Zavala fc Taylor! 12003) with the 
integrated RM removed (Fig. 120k). The 4.6 and 5.1-GHz \ 
values do not fit exactly on the A 2 line, possibly indicating 
a change in the sign of the RM at frequencies lower that 4.6 
GHz. Upstream of the core, the RM map shows a small re- 
gion of very high RM. We were unable to prevent the AIPS 
task RM from unnecessarily rotating the 4.6 and 5.1 GHz 
EVPAs by 180° in this region, hence producing a spurious, 
large RM that is not consistent with the higher frequency 
observations. 

Probing further down the jet in the 7.9-12.9 GHz fre- 
quency range, the sign of the core RM changes to +240 ±90 
rad m -2 (Fig. 120b ). Because we only have three \ values 
in this range, the error is large, but the increase in the 
magnitude of the core RM is significant, and is consistent 
with the general trend of the RM magnitude increasing at 



higher frequencies. Downstream of the core, a band of RM 
of -328 ± 43 rad m" 2 is detected. These bands of RM 
with sign c hanges are a common feature of pc-scale RM 
maps (e.g., IZavala fc Taylor! l200ll ) and are possibly associ- 
ated with bends in the inner jet leading to changes in the 
dominant LoS component of the magnetic field (see Section 
4.5). The jet RM of 162 + 49 rad m -2 again supports a lower 
electron density in the jet compared to the core region, but 
is higher than the jet RM from the lower resolution 4.6-12.9 
GHz data. 

The highest frequency range (15.4-43 GHz) core RM 
of —1008 ± 43 rad m~ 2 changes sign again and increases in 
magnitude (Fig, 120b ). giving a value of a — 1.40 + 0.18. The 
RM fits in the jet are rather poor due to the weak polarized 
jet emission at 22 and 43 GHz. However, there is evidence 
for a significant increase in magnitude of the RM where the 
jet bends. 

The jet EVPAs are aligned with the jet direction after 
removing the effect of Faraday rotation, indicating a trans- 
verse magnetic field structure. There is significant interknot 
polarization and the intrinsic EVPAs remain well aligned 
with the jet direction even as it bends (Fig. 121b ,. b, c). The 
intrinsic core polarization orientation of ~ 40° appears to 
bear no relation to the observed inner jet direction in the 
15-43 GHz xo map, but the higher resolution of our 43- 
GHz map shows a jet component in structural position an- 
gle ~ 210 - 220° from the core (Fig. fl9f). Therefore, the 
EVPAs are also aligned with the jet direction in this region. 

The core degree of polarization decreases from ~ 2.1% 
at 4.6 GHz to a minimum of ~ 1.3% at 7.9 GHz before 
increasing to ~ 3.7% at 22 GHz, with a slight decrease to 
~ 3.5% at 43 GHz. The drop in core polarized flux at 43 GHz 
can be attributed to the newly resolved polarized component 
upstream of the core (Fig. I19f ). The degree of polarization 
in the jet remains relatively constant across the individual 
frequency ranges with m ~ 15-20% at lower frequencies and 
increasing to values of m ~ 25-30% at higher frequencies. 
The degree of polarization maps show an unusual structure 
with the edges of the inner jet having lower m (~ 15%) than 
the central region with m ~ 20% (Fig. [22}. Note that this 
is opposite to the trend expected for a helical jet magnetic 
field, which should yield increased m at the jet edges. This 
unexpected behaviour may be due to turbulence at the edges 
of the jet that disorders the magnetic field. 

There is a substantial increase in m where the outer 
jet bends. Even though the jet direction clos er to the core 
has been observed to vary quite substantially (|Stirlindl2003l ; 
iMutel fc Dennl2005h . the outer jet has always been observed 
to bend to the south-east (e.g., Fig. [22} and the higher m at 
this bend possibly indicates that it is caused by a constant 
interaction with the surrounding medium. 



4 DISCUSSION 

4.1 Core Degree of Polarization 

The core degree of polarization has a V-shaped distribution 
for four (0954+658, 1156+295, 2007+777 and 2200+420) 
of the six sources in this sample, where it initially de- 
creases rapidly from its value at the highest frequency and 
then, after reaching a minimum, increases again toward the 
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Figure 16. Total intensity (/) and polarization (p) maps for 2007+777: (a) 4.6 GHz; / peak 0.55 Jy bm ; lowest contour 1.10 mjy 
bm _1 ; p vectors clipped at 1.0 mjy bm -1 , (b) 7.9 GHz; / peak 0.60 Jy bm -1 ; lowest contour 1.20 mjy bm -1 ; p vectors clipped at 1.0 
mjy bm -1 , (c) 12.9 GHz; / peak 0.54 Jy bm -1 ; lowest contour 1.79 mjy bm -1 ; p vectors clipped at 1.5 mjy bm -1 , (d) 15.4 GHz; I 
peak 0.47 Jy bm -1 ; lowest contour 2.12 mjy bm -1 ; p vectors clipped at 1.5 mjy bm -1 , (e) 22 GHz; J peak 0.38 Jy bm -1 ; bottom 
contour 3.46 mjy bm -1 ; p vectors clipped at 5.5 mjy bm -1 , (f) 43 GHz; / peak 0.39 Jy bm -1 ; lowest contour 6.62 mjy bm -1 ; p vectors 
clipped at 7.5 mjy bm -1 . The / contours increase by factors of two in all maps. Length of p vectors represent relative polarized intensity. 



lowest frequency (Fig. [23}. We investigate whether Fara- 
day depolarization across the observing beam could cause 
the observed de c rease in m using the mo del described in 
I Zavala fe Tavlorl l|2004h and IJorstadl (I2007T ). where m(%) = 
TiRM=o|sinc(RMA 2 )|. We assume that the core degree of 
polarization at 43 GHz is ~ 90% of dirm^o as an initial 
estimate. For 0954+658 and 1156+295, the model has rel- 
atively good success in predicting the first null value in 
m. However, the degree of polarization at 22 GHz falls off 
faster than expected in both cases (see Figs. 123b . b). For 
2007+777 and 2200+420, the model describes the initial 
fall-off much better but underestimates m at 7.9 and 8.9 
GHz, although these values may be affected by blending 



with inner jet components (see Figs. [23b. d). Therefore, we 
conclude that, while Faraday beam depolarization may have 
a significant effect, it cannot completely explain the initial 
fall-off in m. The increase in the degree of polarization at 
low frequencies is most likely due to higher polarized inner 
jet regions becoming bl e nded with the core, a s desc ribed in 
iGabuzda fc Cawthomel l| 19961 ) and lGabuzdal l|l999h . 

It's also worth noting that this distribution is simi- 
lar to the fracti o nal li near polarization spectrum derived in 
IJones fc O'Dell (fl977h . with a turnover frequency (u m ) for 
the core component dominating the observed polarization 
of ~10-15 GHz. In their model the degree of polarization 
is high at frequencies much greater than v m , then decreases 
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Figure 17. Source: 2007+777. (a) 4.6-8.9 GHz RM map (colour scale RM range: -0.5-1 xlO 3 rad m -2 ) with 4.6 GHz total intensity 
contours (/ peak 0.53 Jy bm _1 ; lowest contour 1.33 mjy bm _1 ). The solid arrow points to the core region where the \ vs. A 2 fits shown 
in the inset were found; the degree of polarization is plotted underneath for each corresponding polarization angle. The dashed arrow 
points to the indicated slice taken across the RM map; the plotted black line shows the RM fit at every pixel across the slice with the 
grey area indicating the error in the fit for each point, (b) 12.9-22 GHz RM map (colour scale RM range: -0.5— 1x10 s rad m — 2 ) with 
12.9 GHz total intensity contours (/ peak 0.50 Jy bm" ; lowest contour 2.00 mjy bm" 1 ). All I contours increase by factors of two. Full 
resolution image in journal version. 



rapidly to a minimum just below v m before increasing again 
toward low frequencies, similar to the sort of distribution we 
observe. To properly analyze this effect we need to find the 
spectral turnovers of individual model fitted components in 
the core; this will be investigated in a future paper. With 
regard to the two sources which did not show this V-shaped 
distribution, the core m for 1749+096 may be affected by 
the birth of a new strongly polarized component and for 
1418+546 polarization was not detected at 22 or 43 GHz. 



4.2 Core Rotation Measure 

In all cases, the magnitude of the RM increases at higher fre- 
quencies, consistent with an increase in the electron density 
and/or magnetic field strength closer to the central engine. 
This is quantified using the relation |RM cor c,i/ |oc v a \ see Ta- 
ble [5] for values of a. Table [5] shows that two of our sources 
(1156+295 and 2007+777) are consistent with the value of 
a ~ 2 expected for a stratified jet with a Faraday rotating 
boundary layer sur rounding a con ically expanding jet, simi- 
lar to the results of lJorstadl (|2007f l. However, 0954+658 and 
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Figure 18. Total intensity (/) contours with intrinsic polarization orientation (xo) f° r 2007+777 predicted from (a) 4.6-8.9 GHz RM 
map and (b) 12.9-22 GHz RM map. Same / contour levels as in Fig. 1171 



1418+546 have significantly higher values of a = 3.3 and 3.8, 
indicating much faster electron density fall-offs in the Fara- 
day rotating medium with distance from the central engine. 
In the case of 0954+658, the fit is poor due to a large jump in 
the RM value at high frequencies (using only the two lower 
frequency core RMs we get an estimate of a ~ 1.8), indicat- 
ing that we are observing a region with much higher electron 
density. It is tempting to attribute the jump in core RM from 
15-43 GHz to th e affect of the broad-line region, however, 
IXu fc Cad | 20061 ) derive an upper limit on the size of the 
broad-line region for 0954+658 of ~ 0.4 pc which is much 
smaller than the region probed by our 15-43 GHz observa- 
tions of ~ 2 pc. Another possibility is that the steep increase 
in the core RM is relat ed to steep press ure gradients in the 
region confining the iet. lLobanovl ljl998t ) showed that the jet 
particle distribution follows approximately the same power- 
law as the external pressure at distances much greater than 
where the j et becomes supersonic using t he hy drodynamic 
jet model of lGeorganopoulos fc Marscherl (|l998l ). Hence, the 
large values of a may correspond to jets confined by ther- 
mal pressure and/or ram pressure from a spherical disk- 
wind outflow surrounding a conically expanding jet (e.g., 
iBogovalov fc Tsinganosl [2005), while lower values of a (e.g., 
a = 1.4 for 2200+420) could be related to highly collimated 
outflows possibly with a signi ficant contribution to c onfin- 
ment from magnetic pressure l|Komissarov et al"1l2007l ). 



The intrinsic core RM appears to be anti-correlated 
with the core degree of polarization (Fig. I24[). with a dearth 
of high polarized cores with large intrinsic RMs (with the 
exception being 2007+777 in the high frequency range where 
we observe an optical depth transition for the core polariza- 
tion). This is consistent with the findings of lZavala fc Tavlorl 
(|2003t ) from a much larger sample for both quasars and BL 
Lac objects. They attribute this to depolarization from a 
gradient in Faraday rotation across the beam, but as we 
have seen this alone cannot completely explain the observed 
frequency dependence of m. 



4.3 Intrinsic Polarization Orientation 

After correction for Faraday rotation, the jet polarization 
orientation is approximately aligned with the jet direction 
for 0954+658, 2007+777 and 2200+420, and approximately 
transverse to the jet direction for 1418+546, while 1156+295 
displays a spine-sheath polarization structure. EVP As that 
are aligned with the jet direction have often been interpreted 
as being associated with transverse shocks that enhance 
the magnetic field component in the plane of compression 
(e.g., iHughes. Aller fc Allerl Il989h . but this interpretation 
has difficulty explaining significant aligned inter-knot polar- 
ization and aligned EVPAs fol l owing jet bends (e.g., Fig. l21[) . 
iLvutikov. Pariev fc Gabuzdal l|2005t ) have shown that this 
bi-modal distribution of the jet EVPAs is consistent with 
the presence of a helical magnetic field geometry, where the 
aligned EVPAs correspond to the dominant toroidal com- 
ponent of a relatively tightly wound helical field and the 
transverse jet EVPAs correspond to the poloidal component 
of a relatively loosely wound helical magnetic field. Regions 
of longitudinal magnetic field near the jet edges (Fig. [TJi) 
are often attributed to int eraction between the jet and the 
surrounding medium fe.g.. lAttridge et alj|2005h . but in fact 
can also be explained as the projected longitudinal compo- 
nent of a helical magnetic field at the jet edges, when there 
is sufficient resolution across the jet. 

The intrinsic core EVPAs are more complicated to anal- 
yse due to possible blending of polarized components in 
the innermost jet at low frequencies. They also vary with 
frequency-interval because the different core RMs can im- 
ply different zero-wavelength angles, due to the different 
regions being probed. 2200+420 is the only source whose 
intrinsic core EVPA maintains an approximately constant 
direction, with xo(4.6-12.9GHz) = 37 + 4° and Xo(15.4- 
43GHz) =46+1°, which is aligned with the inner jet direc- 
tion seen at 43 GHz implying a dominant transverse mag- 
netic field all along the pc-scale jet. For all other sources, 
we use the highest frequency-interval core \o value or, if 
the RM fits are unreliable or the core region is affected by 
blending, we use the 43-GHz core polarization angle. For 
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Figure 19. Total intensity (7) and polarization (p) maps for 2200+420: (a) 4.6 GHz; / peak 1.43 Jy brn" 1 ; lowest contour 3.58 mjy 
bm _1 ; p vectors clipped at 1.5 mjy bm" 1 , (b) 7.9 GHz; / peak 1.83 Jy bm _1 ; lowest contour 1.83 mjy bm" 1 ; p vectors clipped at 1.5 
mjy bm -1 , (c) 12.9 GHz; J peak 1.90 Jy bm -1 ; lowest contour 3.22 mjy bm -1 ; p vectors clipped at 1.8 mjy bm -1 , (d) 15.4 GHz; I 
peak 1.78 Jy bm -1 ; lowest contour 3.39 mjy bm -1 ; p vectors clipped at 2.0 mjy bm _I , (e) 22 GHz; J peak 1.73 Jy bm -1 ; bottom 
contour 6.06 mjy bm -1 ; p vectors clipped at 4.0 mjy bm -1 , (f) 43 GHz; I peak 1.90 Jy bm -1 ; lowest contour 9.14 mjy bm -1 ; p vectors 
clipped at 5.5 mjy bm -1 . The / contours increase by factors of two in all maps. Length of p vectors represent relative polarized intensity. 
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Figure 20. Source: 2200+420. (a) 4.6-12.9 GHz RM map (colour scale RM range: -0.5-0.7xl0 3 rad m~ 2 ) with 4.6 GHz total intensity 
contours (/ peak 1.39 Jy bm _1 ; lowest contour 4.18 mjy bm _1 ). The solid arrow points to the core region where the \ vs - ^ 2 fits shown 
in the inset were found; the degree of polarization is plotted underneath for each corresponding polarization angle. The dashed arrow 
points to the indicated slice taken across the RM map; the plotted black line shows the RM fit at every pixel across the slice with the 
grey area indicating the error in the fit for each point, (b) 7.9-12.9 GHz RM map (colour scale RM range: -0.7-0.6xl0 3 rad m~ 2 ) with 
7.9 GHz total intensity contours (/ peak 1.80 Jy bm" 1 ; lowest contour 3.60 mjy bm" 1 ). (c) 15.4—43 GHz RM map (colour scale RM 
range: -2.5-1.5xl0 3 rad m -2 ) with 15.4 GHz total intensity contours (I peak 1.78 Jy bm -1 ; lowest contour 3.56 mjy bm -1 ). All / 
contours increase by factors of two. Full resolution image in journal version. 
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Figure 21. Total intensity (/) contours with intrinsic polarization orientation (xo) for 2200+420 predicted from (a) 4.6-12.9 GHz RM 
map, (b) 7.9-12.9 GHz RM map and (c) 15.4-43 GHz RM map. Same I contour levels as in Fig. [20] 



0954+658, X43GHz ~ 30° which is approximately transverse 
to the inner jet direction of ~ —50° . It is unclear whether the 
magnetic field is parallel or perpendicular to the core EVPA 
since we are not sure if the polarized emission is coming pre- 
dominately from an optically thick or thin emitting region 
within the observed VLBI "core". In the case of 1156+295, 
the value of Xo(12.9 - 43GHz) = -17 ± 4° is offset from 
the observed inner jet direction of ~ 0°, but a jet direction 
of ~ —17° at higher frequencies would n ot be i ncons istent 
with the helical jet trajectory proposed bv lHonj (|2004T ) . The 
highest frequency RM correction for the core of 1418+456 
gives xo(7.9 — 15.4GHz) = 109 ± 3° compared to a inner 
jet direction of ~ 130°. The approximate 90° change in the 
EVPA between the core and inner jet suggets that the ob- 
served core polarization at high frequencies may be from the 
optically thick region, whereas the observed core polariza- 
tion at lower frequencies is dominated by polarized emission 
from the optically thin inner jet region. In this case, the 
magnetic field direction would be approximately longitudi- 
nal in both the core and jet. The core and inner jet spec- 
tral indices for 1418+546 of +0.3 and —0.6, respectively, 
support this hypothesis, as well as the degree of polariza- 
tion at 15 GHz, which increases from ~ 2% in the core to 
~ 10% in the inner jet. The 43-GHz core polarization angle 
for 1749+096 of ~ —50° is perpendicular to the jet direction 
of ~ 40°. It's unclear whether the dominant magnetic field 
structure is transverse or longitudinal in this source, since 
the EVPAs remain approximately transverse to the jet di- 
rection across the entire frequency range, while the spectral 
index indicates that the core is strongly self-absorbed at low 
frequencies (a ~ +0.8) but turns over at high frequencies 
with a ~ —0.1. Using the 43-GHz core polarization angle 
for 2007+777, we find the core EVPA of ~ -80° approxi- 



mately aligned with the jet direction of —90°. This indicates 
that this source is dominated by a transverse magnetic field 
structure in both its pc-scale core and jet. 

4.4 Transverse Jet Structure 

Analysis of the transverse RM structure of 0954+658, 
1418+546 and 1156+295 reveals the presence of positive 
asymmetric gradients in RM across the jet, which are strong 
signatures for the presence of a helical magnetic field ge- 
ometry in the Faraday rotating medium. RM gradients in 
0954+658 and 1156+ 2 95 ha ve pr eviously been reported by 
iMahmud fc Gabuzda! ((2008) and I Gabuzda^ et al l l|2008h on 
similar scales, and 0954+658 also has an RM gradient on 
decaparsec scales (jHallahan fc Gabuzdal [20081 ) . For a per- 
fectly side-on view of a helical magnetic field, one would 
expect a symmetric gradient, with the LoS magnetic field 
changing sign acr oss the jet. However, as was shown in 
lAsada et alj (|2002l ). a positive asymmetric gradient (as de- 
tected in 0954+658, 1418+546 and 1156+295) can be under- 
stood if the viewing angle in the jet rest frame is less than 
the pitch angle of the helix. All the observations show that 
the RM gradients in these three sources maintain a constant 
direction across the jet. 

An excellent diagnostic of the internal magnetic field 
structure are transverse profiles of the synchrotron emission 
from the jet. However, determining the "ridge line" along the 
jet from which to take orthogonal slices can prove rather dif- 
ficult. Using model fitted total intensity components along 
the jet is not very useful because some jets have quite broad 
funnels through which components move following different 
paths. Since the transverse jet emission itself can be asym- 
metric, using the peak intensity with distance along the jet 
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Figure 22. 7.9 GHz total intensity and degree of polarization 
(m), ranging from to 45 percent, for 2200+420 (/ peak 1.82 Jy 
bm" 1 ; bottom contour 2.74 mjy bin -1 ). The sliced region shows 
the symmetric distribution of m across the jet where it has a high 
m along the spine of the jet and lower m at the edges. 



is also not always appropriate. One option is to estimate the 
jet edges and define the ridge line as the mid point. How- 
ever, estimating the jet edges where the emission is weak 
is often difficult and error prone. For the transverse pro- 
files presented in this paper, we defined the normal to jet 
direction "by eye" and then checked that the profiles did 
not change dramatically for small changes of the slice direc- 
tion, so that our conclusions didn't depend on our choice of 
profile. 



Asymmetric distributions of transverse profiles of total 
intensity and polarization appear to be a common feature 
in pc-scale jets (e.g., Fig. I25|l . Helical magnetic fields pro- 
duce asymmetric total intensity and polarization transverse 
profiles unl ess the jet is viewed at exactly 90° in the jet 
rest frame ljLaind[l98l| ). Models of a magnetic field struc- 
ture in a cylindrical jet containing a uniform helical field 
with a tangled component (to reduce the degree of polariza- 
tion to observed levels) are able to reproduce the observed 
offsets between the peaks in the total and polarized inten- 
sities (Papageorgiou & Cawthorne 2008, in prep.) providing 
strong support for a helical magnetic field in the synchrotron 
emitting plasma. 

Furthermore, the total intensity is maximum where the 
component of the magnetic field is perpendicular to the LoS 
(/ tx Sin[#]), while the RM is maximum where the com- 
ponent of the magnetic field is parallel to the LoS (RM oc 
Cos[#]). Hence, for a helical magnetic field viewed at an an- 
gle other than 90° in the jet rest frame, the magnitude of the 
RM should be greatest on the opposite side of the jet from 
the total intensity peak. In the case of 0954+658 (Fig. 1251 
a) and 1156+295 (Fig. 1251 b), it is not clear as the offset 
from the centre of the jet is quite small, ~ 6% and ~ 3% 
of the their respective beam-sizes, but this is not the case 
for 1418+546 (Fig. 1251 c) where the offset from the centre 
of the jet to the I peak is ~ 20% of the beam. Since we do 
not expect the thermal plasma causing the observed Fara- 
day rotation to be mixed in with the synchrotron emitting 
material, the RM and I slice may not correspond to exactly 
the same helical field structure. If the jet is electromagneti- 
cally dominated, then the field lines must close back in the 
accretion disk, and it's possible that the toroidal component 
of the field in the sheath surrounding the jet is orientated 
in the opposite direction, causing the magnitude of the RM 
to be largest on the same side of the jet as the peak offset 
of the synchrotron emission (Figs. 1251 b and c). However, 
comprehensive modelling of these profiles is clearly required 
to reach any firm, quantitative conclusions. 

4.5 Core RM Sign Changes 

RM sign changes in different frequency-intervals are ob- 
served in 0954+658, 1418+546 and 2200+420. In the case of 
0954+658 and 2200+420, the dominant jet magnetic field is 
transverse throughout the jet (even as the jet bends, in the 
case of 2200+420) implying the presence of a global trans- 
verse magnetic field structure. If we suppose that a helical 
magnetic field surrounds the jet, the observed RM sign re- 
versals can be explained by slight bends of the relativistic 
pc-scale jet, for example, due to a collision with material in 
the paren t galaxy or some instabilities inherent in the jet 
itself fsee lAguddl2008h . 

To show this, we generated simple model transverse 
RM profiles considering a helical magnetic field threading 
a sheath of Faraday rotating material surrounding the jet. 
On the left panels of Fig. 1261 the top side of the jet, which 
corresponds to negative values on the transverse profile axis, 
has the toroidal component of the helical field orientated to- 
wards the observer. Relativistic aberration of the radiation 
emitted by the jet means that radiation emitted at 90° in 
the jet rest frame is observed at an angle (ff) of ~ 1/F in the 
observers frame, where V is the bulk Lorentz factor of the 
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Figure 23. Core degree of polarization vs. A 2 from 43-7.9 GHz for (a) 0954+658 and (b) 1156+295, and from 43-4.6 GHz for (c) 
2007+777 and (d) 2200+420. Solid line shows predicted beam depolarization. See Section 4.1 for details. 



jet. Hence, a side-on view (6 ~ 1/r) of a helical magnetic 
field (Fig. 1261 top left) will have an RM that varies uniformly 
across the jet from positive to negative; considering an un- 
resolved transverse jet structure, a zero net RM should be 
observed. If the jet bends towards us, we have a "head-on" 
view of the helical magnetic field (i.e., < 1/r), providing 
an asymmetric RM profile across the jet (Fig. 1261 middle 
left). The dominant LoS magnetic field component comes 
from the top half of the jet and, for a beam-size compara- 
ble to or bigger than the diameter of the jet, the transverse 
RM structure is resolved out and a net positive RM will be 
observed in this case. Conversely, for a "tail-on" view of a 
helical magnetic field (i.e., 8 > 1/r) (Fig. 1261 bottom left), a 
negative RM will be observed. (A longitudinal jet magnetic 
field with a change in the angle to the LoS could also cause 
the observed effects, but we have considered a helical field 
here, because the observed intrinsic EVPAs usually indicate 
a dominant transverse magnetic field and all three sources 
display RM gradients that provide direct evidence for helical 
fields associated with their jets.) 

The RM sign changes could also be due to accelera- 
tion/deceleration of the jet since this will change F and 
therefore, the viewing an gle relative to 1/r. Recent observa- 
tions (e.g- Jjorst ad 2005) have pr ovided evidence for pc-scale 
acceleration in blazar jets, and IVlahakis fe Konlgl (|2004[ ) 
show that an extended (i.e., pc-scale) acceleration region is 
a distinguishing characteristic of MHD driving of relativistic 
outflows compared to purely hydrodynamic models. Impor- 
tantly, our explanation for the RM sign changes implies rela- 



tivistic bulk motions of the Faraday rotating material in the 
sheath, at least in the innermost part of the jet. This would 
support any theoretical simulations indicating a mildly rela- 
tivistic outflow from the accretion disk surrounding the radio 
emitting part of the jet. 

Therefore, it's possible that either pc-scale bends or an 
accelerating/decelerating jet surrounded by a helical mag- 
netic field could be responsible for the observed RM sign 
changes, because both effects cause the viewing angle to 
change relative to the angle 1/r along the jet. As noted 
in Section 2.1, VLBI resolution is usually not sufficient to 
completely resolve the true optically thick core, so on scales 
smaller than the observed VLBI core, "core" RMs with dif- 
ferent signs could be derived from observations at different 
frequencies (i.e., probing different scales of the inner-jet). 



5 CONCLUSIONS 

We have analysed the polarization and RM distribution of 
six blazars from 5-43 GHz and a summary of our results is 
outlined below. 

1. We find that the magnitude of the core RM in- 
creases systematically with increasing frequency and is well 
described in most cases by |RM coro |(x v a . The values of a 
vary from 0.9 to 3.8. Intermediate values of a ~ 2 would be 
expected for a boundary layer of Faraday rotating material 
surrounding a conically expanding jet. Low values of a may 
correspond to highly collimated outflows while large values 
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Figure 24. Plot of intrinsic core RM (RM; ntI .; ns i c = 
RMobacrvod (1 + z) 2 ) for all sources vs. the corresponding core 
degree of polarization (m) listed in Table \5\ There is clearly a 
lack of highly polarized cores with large core RMs. 



of a are possibly related to large external pressure gradients 
from a spherical disk- wind type geometry. 

2. The core degree of polarization for our observed 
sources decreases rapidly from its value at 43 GHz, to a 
minimum somewhere between 8-15 GHz before increasing 
again at lower frequencies. The initial decrease cannot be 
completely explained by Faraday beam depolarization, and 
we attribute the low frequency increase to blending of higher 
polarized inner jet components. 

3. We detect gradients in the RM across the jets of 
0954+658, 1156+295 and 1418+546, which provides evi- 
dence for the presence of helical magnetic fields in a sheath 
or boundary layer surrounding the jet. 

4. After correcting the polarization vectors for the ef- 
fect of Faraday rotation, we find the jet EVPAs to be ap- 
proximately aligned with the jet direction for 0954+658, 
2007+777 and 2200+420. 1156+295 displays a spine-sheath 
polarization structure, and in 1418+546 the intrinsic jet EV- 
PAs are approximately transverse to the jet direction. The 
EVPAs of 2200+420 display a continuous structure that re- 
mains aligned with the jet direction even as it bends. This 
shows that the magnetic field structure in the synchrotron 
emitting plasma is dominated by an ordered transverse com- 
ponent in all of these sources except for 1418+546. A heli- 
cal magnetic field geometry can neatly explain both the bi- 
model distribution of the jet EVPAs and the well ordered 
polarization structures. 

5. For three of the five sources with detected Faraday 
rotation (0954+658, 1418+546 and 2200+420), we find that 
the core RM also changes sign with distance from the central 
engine. We provide an explanation for this by considering a 
boundary layer of Faraday rotating material threaded by 
a helical magnetic field, with bends in the relativistic jet or 
acceleration/deceleration of the jet leading to changes in the 
viewing angle relative to the angle 1/F (corresponding to a 
viewing angle of 90° in the jet rest frame). This in turn gives 
rise to changes in the sign of the dominant LoS component 
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Figure 25. (a) Transverse profiles of the 4.6 GHz total intensity 
and polarized intensity of 0954+658 along with the 4.6-15.4 GHz 
RM slice across the region indicated in Fig. [5^. (b) Transverse 
profiles of the 4.6 GHz total intensity and polarized intensity of 
1156+295 along with the 4.6-15.4 GHz RM slice across the region 
indicated in Fig. |8^.. (c) Transverse profiles of the 4.6 GHz total 
intensity and polarized intensity of 1418+546 along with the 4.6- 
8.9 GHz RM slice across the region indicated in Fig. Illfa . 
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Figure 26. Top right: Cartoon of a conical jet with the heavy line indicating a transverse profile of the jet. The circle is an example of 
the finite resolution of our VLBI observations. Middle right: Due to the relativistic motion of the jet towards us, radiation emitted at 
< 90° in the jet rest frame is observed at < 1/r where T is the bulk Lorentz factor of the jet. Bottom right: Radiation emitted at > 90° 
in the jet rest frame is observed at > 1/r. Top left: Transverse RM profile of a jet surrounded by a sheath of Faraday rotating material 
and viewed at 90° in the jet rest frame. Due to the finite resolution of our observations a zero net RM would be observed. Middle left: 
Transverse RM profile for a jet viewed "head-on" (< 90° in the jet rest frame). If we cannot resolve the transverse jet structure a net 
positive RM will be observed. Bottom left: Transverse RM profile for a jet viewed "tail-on" (> 90° in the jet rest frame). In this case, a 
net negative RM will be observed if the transverse structure is resolved out. 



of the magnetic field, and thereby to changes in the sign of 
the RM. 
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